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ABSTRACT 

We suggest that the unique X-ray source 1E161348-5055 at the centre of the 
supernova remnant RCW 103 consists of a neutron star in close orbit with a 
low mass main sequence star. The time signature of 6.67 hr is interpreted as the 
neutron star's spin period. This requires the neutron star to be endowed with 
a high surface magnetic field of ~ 10 15 G. Magnetic or/and material (propeller) 
torques are able to spin rapidly the young neutron star down to an asymptotic, 
equilibrium spin period in close synchronism with the orbital period, similarly to 
what happens in the Polar Cataclysmic Variables. 1E161348-5055 could be the 
first case of a magnetar born in a young low mass binary system. 

Subject headings: Stars: binaries: general - Stars: magnetic fields - Stars: 
neutron - Stars: cataclysmic variables - Stars: supernovae: individual: RCW 103 
- X-rays: binaries 



Introduction 



The soft X-ray source IE 161348-5055 (hereafter IE) was discovered by lTuohy and Garmire 



( 119801 ) with the Einstein observatory, close to the geometrical centre of the young supernova 
remnant (SNR) RCW 103. Owing to its thermal-like spectrum and to the lack of a counter- 
part at radio or optical wavelengths, IE was classified as the fi rst example of a ra dio-quiet, 
isolated neutron star. The host S NR is very young (~ 2000 vr. ICarter et al.lll997l ) and it is 
located at a distance of ~ 3.3 kpc ([Caswell et al.lll975l ; iReynoso et al.ll2004f ) . The association 
of IE and RCW 103 is very robust. The point source lies within 1 5" of the SNR centr e; 
moreover, IE and RCW 103 have consistent distance measurements ( IReynoso et al.l 120041 ). 
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The original interpretation of IE as an isolated neutron star was later questi oned by the 
observation of a large variability on few years' time scale (jGotthelf et al.lll999l ). This was 
confirmed in recent years, when a factor ~ 100 brightening (from ~ 9 x 10~ 13 erg cm" 2 s" 1 
to ~ 7 x 10 -10 erg cm~ 2 s _1 ) was discovered, comparing the first and the second Chandra 
observations of the so urce, performed in September 1999 and in March 2000, respectively 
(IGarmire et al.ll2000bl ). Even more puz zling, a possible long periodicity at ~ 6 hours was 
hinted by the first Chandra observation (IGarmire et al.ll2000al ). but it w as not recognised in 
the second Chandr a data set. Subsequent observ ations with Chandra (jSanwal et al.ll2002l ) 
and XMM-Newton (IBecker and Aschenbachl |2002| ) failed to ultimately settle the periodicity 



issue. 



A breakthr ough came with the deep (90 ks) observation of IE performed by XMM- 



Newton in 2005. iDe Luca et al.l (120061 ) caught the source in a low state (1.7x 10 12 erg cm 2 s 



and reported conclusive evidence for a strong (~ 50%), nearly sinusoidal modulation of IE 
at P = 6.67 ± 0.03 hr. The source spectrum, thermal-like and well described by the sum 
of two black bodies, varies along the 6.67 hr cycle and is harder at the pulse maximum 
( IDe Luca et al.l 120061 ). The 6.67 hr periodicity was also recognised in the former XMM- 
Newton data set (collected in 2001, when the source was a factor ~ 5 brighter), although 
with a much smaller pulsed fraction (~ 12%) and with a remarkably different, complex light 
curve, featuring two narrow minima (or "dips") per period, separated by 0.5 in phase. No 
faster periodicities were found, with a 3 a upper limit of 10% on the pulsed fraction down 
to P = 12 ms. 



De Luca et al.l (120061 ) reported also a complete picture of the peculiar long-term vari- 



ability of IE. The source has been continuously fading along more than 5 years, since the 
1999-2000 re-brightening. In August 2005 IE was seen with a flux similar to the pre-outburst 
level. XMM-Newton data clearly show that the spectrum of the source is also evolving. 
While a double black body model yields in any case the best description of the spectrum, 
the emission is harder and more absorbed when the source is brighter. 

On the optical side, deep observations of the field with the ESO/VLT and with HST 
showed nothing b ut two /three very fai nt IR sources (H ~ 22) consistent with the accurate 
Chandra position (IDe Luca et al.ll2006l . and reference therein). No firm conclusions could be 
drawn about the possible association of any of such sources with IE. 
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1.1. Isolated vs Double Star Scenario 

The combination of long-term variability, 6.67 hours periodicity, young age and under- 
luminous optical/IR counterpart makes IE a unique source among compact objects in SNRs. 



As discussed by iDe Luca et al.l (120061 ). IE could either be a very young low mass X-ray 
binary (LMXB) system, featuring a 2000 yr old compact object and a low mass companion 
star (M4, or later) in an eccentric orbit, or a very peculiar isolated neutron star. In the 
latter case, IE could only be a magnetar, dramatically slowed down, possibly by a propeller 
interaction with a debris disc. Both scenarios require highly non-standard assumptions. 

The binary scenario is appealing since it provides an easy explanation of the long 6.67 hr 
periodicity, as due to the orbital period of the system. However, the observed properties of IE 
are dramatically different from those of any known LMXB. The system is orders of magnitude 
dimmer than persistent LMXBs. Moreover, a very peculiar "double accretion" mechanism is 
required to explain the unusual, large flux and spectral variations along the 6.67 hr cycle as 
an orbital modulation, as well as the dramatic long term variability. Indeed, wind accretion 
along an eccentric orbit could naturally explain the quasi-sinusoidal modulation observed in 
the low state. However, the light curve of the source in its "active state" is very complex. 
In order to explain its features within such a scenario - as due to varying occultation of the 
central source - the existence of very peculiar structures at the rim of a (possibly transient) 
accretion disc is required. 

The isolated magnetar scenario could easily explain the phenomenology of the periodic 
modulation, with the star rotation bringing into view or hiding different emission regions. 
Indeed, the light curve profile is very similar to the one observed in most magnetars. The 
long-term variability (both in luminosity and in pulse profile) is also reminiscent of the one 
seen in a few "transient" anomalous X-ray pulsars (AXPs). However, all known magnetars 
spin thousands of time faster than IE, with periods well cluste red in the 5 — 12 s r ange. 



The spin history of IE should have been dramatically different. IDe Luca et al.l ((20061) and 



Lil (120071 ) showed that propeller interaction with a supernova debris disc could provide a 
very efficient slowing-down mechanism. Such a mechanism could brake the neutron star to 
6.67 hr in 2000 yr, provided that the star was born with a ~ 10 15 G magnetic field and an 
initial period not faster than 300 ms, in order to avoid an initial "eject or" phase which would 



have pushed away any surrounding material (see e.g. lLipunovlll992l ). Indeed, evidence for 



a deb ris disc surrounding the AXP 4U 0142 + 61 has been recently obtained (jWang et al. 



20061 ). However, this source spins at P = 8.4 s, a very typical value for an AXP. Thus, 
peculiar birth properties, coupled to some fine tuning, should be required to explain the 
unique spin evolution of IE in the isolated magnetar scenario. 
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In this paper we will explore an alternative scenario, also suggested by iPopovl (120061 ) 
as a possibility, for the nature of IE. We assume that IE is a binary system, hosting a 
recently born neutron star and a low mass companion. Within such a frame, we will suppose 
the 6.67 hr periodicity to be the spin period of the neutron star. The required spin-down 
mechanism is provided by material and/or magnetic interaction between the neutron star 
and its stellar companion. A magnetar-like field of ~ 10 15 G is required for these interactions 
to brake the neutron star's rotation that slow. 

The binary system envisaged in our model has some interesting analogies with some 
magnetic Cataclysmic Variables. Inde ed, the present i nvest igation is partially motivated by 
the similarity, already pointed out by iDe Luca et al.l (120061 ). of several temporal properties 
of IE (the period value, the pulse shape and its long-term variability) with those of Polars 
(also k nown as AM H erculis), a class of strongly magnetised Cataclysmic Variables (CVs; 
see e.g. ICropperlll990l for a review). 



In such systems the magnetic interaction between the components is so strong to lock the 
white dwarf in synchronous rotation with the orbit, with periods up to few hours (the longest 
period of a Polar, V1309 Ori, is of 8 hr, iGarnavich et al.lll994j ). In addition, the magnetic 
pressure prevents the formation of a disc around the white dwarf, and the accretion flow 
is directly channelled from the inner Lagrangian point Li to one or both the white dwarf's 
poles. Our model owes a good deal of inspiration to the current theoretical models of strongly 
magnetic CVs. The reader is warned, however, that the "Polar" model we envisage for IE 
may differ in some important details from its CV counterpart. 

We shall address the fundamentals of the model, deferring to a further paper the detailed 
explanation of the light curve and the spectral properties of the source. 

In Section [2] we elaborate more quantitatively on the "Polar" model. In Section [3] 
we discuss the dynamical equations which rule the system's evolution after the supernova 
explosion which left the SNR RCW 103. In Section H] we solve numerically these equations, 
and briefly discuss the properties of their asymptotic equilibrium solutions (this topic is 
covered more fully in Appendix [AJ. In Section [5] we cover a number of important issues 
such as the secondary's evaporation and the future evolution of IE. Our conclusions are 
summarised in Section [HI 



2. Setting up the Model 



We assume that IE is a binary system made by a neutron star (NS) of mass Mi = 1.4 M® 
and a low main sequence star (a red dwarf, RD) of mass M2 = 0.4 M® (consistent with the 
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current optical/infrared constraints). This system survived the supernova explosion which 
led to the birth of RCW 103. The orbit is expecte d to be born quite eccentric, with e 
most likely in the range 0.2 — 0.5 (IDe Luca et al.ll2006l ); this eccentricity deserves some com- 
ments. In a narrow system like that we envisage for IE, the tidal torques are expected 
to be quite effective in synchronising the companion star's rotation with the orbital period 



and in circularising the orbit. Unfort unate 
ciency of these two proces ses (see e.g. IZahn 



P 
il 



y, there is no general consensus about the effi 



1977 



Tassoujll988l : lRieutord and Zahnl fll997h : 



Tassoul and Tassoull 119971 : see also iMeibom and Mathieul 120051 and references therein) . The 
most conse rvativ e estim ates of the synchronisation and circularisation time scales are those 
derived bv lZahnl Jl977h (his Equations (4.12) and (4.13)) 
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where L 2 ~ 2 x 10 32 erg s _1 is the main sequence star's luminosity, q = Mi/M 2 is the mass 
ratio, I 2 and ho. are the m ain sequence star's momentum of inertia and the apsidal motion 



Zahn 



9771 for mo re details). The quantities k 2 and I 2 /M 2 R 2 are both of the 



Motall952l ). and therefore should approximately cancel each other, so a 

10 2 yr 



6.67 hr, M l = 1.4 M and M 2 = 0.4 M would take t 



sync 



constant (see 
same order (see e.g. 
stellar system with P or b 
to synchronise and t C j rc ~2x 10 3 yr to circularise. The companion star is thus expected to 
be synchronous, but the orbit may still be eccentric. 

Although the eccentricity may contribute to modulate the X-ray pulsations of IE, its 
introduction at the present level of discussion would complicate the model without any real 
benefit. As it will be discussed in Sectional our results are not expected to be too sensitive 
to the use of a circular orbit approximation. For this reason, we suppose that the orbit is 
circular, with the two stars separated by the fixed distance a. If the NS's period Pi = 6.67 hr 
is equal (or close to) the orbital period P or b, a is given by Kepler's third law 

2/3 / M \ 1/3 



1.5 x 10 



11 



P, 



orb 



6.67 hr 

where M = Mi + M 2 is the binary's total mass. 



1.8 M 



cm, 



(3) 



We assume that the interaction between the stars and the accretion flow are dominated 
by the star's magnetic field, as happens in canonical Polar CVs. This requires the condition 
Ra > a, where 

-2/7 

1.8 x 10- (-J^-Y ' l-^-) I - : en, (4) 



10 33 G cm 3 



1Q13 gs-l 



1.4 Mr. 
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is the NS Alfven radius. The mass accretion rate Mi has been normalised to the low state 
X-ray luminosity Lx of IE according to the well-known formula Lx = G Mi Mi/Ri, being 
Ri = 10 6 cm the NS radius. It is easy to check that the requirement Ra > a ~ 10 11 cm 
implies a NS magnetic moment /ii ~ 10 33 — 10 34 G cm 3 , or a magnetar-like surface magnetic 
field Bi ~ fii/Rf ~ 10 15 - 10 16 G. 

The Roche lobe of the NS is entirely filled by the star's own magnetosphere, and under 
this circumstance the formation of an accretion disc is impossible: any mass flow from the RD 
through Li is immediately channelled along the NS's field lines and accreted. Independently 
of the accretion pattern, the balance between the material accretion spin up torque and the 
magnetic spin down propeller torque fixes the equilibrium condition R co ~ Ra between the 
Alfven radius and the co-rotation radius 

R co = (GMi/uf) l,Z , (5) 

where U\ = 2n/Pi is the NS rotation angular frequency. The rotation equilibrium condition 
R co ~ Ra sets the NS final spin period 

Pl -> P - " 6 ( lQ33Gcm3 ) 6/7 ( axJ^'gs-i ) {l^k) (6) 



In a Polar system (and also in IE) the RD is endowed with a large magnetic moment. 
This may be either induced by the compact star or intrinsic. This last alternative is quite 
possible, as it i s known that M dw arf stars may harbour sizeable magnetic fields, even of 
few kG (see e.g. iDonati et al.ll2006l ). and magnetic moments above 10 34 G cm 3 . In this paper 
we assume for the RD the magnetic moment /i2 = 10 34 G cm 3 . 

The X-ray emission from Polars is powered by mass accretion from the secondary. In 
our case, this may still be true, but some important contributions may come from other 
sources as well. Currently known magnetar candidates (Anomalous X-ra y Pulsars, AXPs 
and Soft Gamma-ray Repeaters, SGRs: see IWoods and Thompson! 120061 for a review) are 
observed to emit X-ray radiation with luminosities ranging from ~ 10 33 to ~ 10 36 erg s _1 . 
This emission is believed to be powered by the decay of their ultra-strong magnetic field; 
this mechanism might be also responsible for the ~ 10 33 — 10 35 erg s -1 X-ray luminosity of 
IE. In addition, the spectrum and the long term variability of IE are remarkably similar 
to the characteristics of so me AXPs and SGRs (for example, the AXP XTE J1810-197, 
Gotthelf and Halpernll2005l ). Therefore, the X-ray emission of IE might be (completely or 
partially) explained by typical magnetar emission rather than by accretion, while the main 
role of the companion star invoked in this paper is to explain the discrepancy between the 
6.67 hr pulsation period of IE and the much shorter periods 5 — 12 s observed in AXPs and 
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SGRs. If this is the case, the X-ray luminosity Lx only provides an upper limit to the mass 
accretion rate M\ on the NS: in the low state of IE, for instance, M\ < 10 13 g s _1 . Of course, 
the companion star may transfer mass to the compact object only provided that its size is 
of the sa me size as its Roche lobe. Accor ding to the similarity relation valid for low ZAMS 
stars (e.g. iKippenhahn and Weigertlll994l ) the radius of the companion star is 



R-2 



3.3 x 10 



10 



Mo 



0.4M R 



0.8 



cm, 



(7) 



and is compara ble to the star's Roche lobe radius Rl ~ 4 x 10 10 cm, as computed from 
Equation (2) of lEggletonl (119831 ) . Since R 2 is close to Rl, in view of the uncertainties, we 
need to explore both the possibilities that the system is detached (R 2 < Rl) or that there is 
a Roche lobe overflow from the RD to the NS (R 2 > -Rl). The geometry of the system (i.e. 
detached or not), plays a fundamental role in the NS spin evolution. Indeed, Polars may 
reach exact synchronism provided that th e magnetic torque is stronger than the accretion 
spin-up torque oc Mi (e.g. ICampbelllll986l ) . If most of the X-ray luminosity of IE is powered 
by magnetar emission and not by accretion, Mi is small (the binary may be detached), so 
the magnetic torques are dominant and may synchronise the system. Note that in order to 
attain synchronism the magnetic torque must not (only) be stronger the accretion torque 
now, but it must have been so also over a significant fraction if the system's life. Incidentally, 
we note that if deeper observations will show the mass of the companion to be significantly 
smaller than its current upper limit (0.4 M ), then the companion is most likely to lie well 
within its own Roche lobe. In this case the system is detached, and the X-ray luminosity is 
clearly due to the magnetar emisson. On the other hand, if the accretion torque is strong 
enough, it may keep the NS away from exact synchronism. In this case IE should be more 
similar to an Intermediate Polar (IP). The spin evolution of these systems is dominated by 
the material torques, and their action leads to rotation equilibri um periods distribu ted over 
a fairly wide range of values, all shorter than the orbital period (jNorton et al.l 12004 ). 



In the next Section we set up all the formal machinery to deal with both these cases, 
aiming to follow the spin evolution of IE. 



3. Dynamics of the System 

In this Section we focus on the evolution of the NS spin and the orbital period after the 
supernova explosion. In order to keep our analysis as simple as possible, we assume that all 
the torques are directed perpendicular to the orbital plane. Any misalignment would lead 
to a rotation of the orbital angular momentum: this just complicates the equations, without 
any relevant benefit for the physical insight into the problem. 
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We evaluate the torques acting on the NS, as well as those acting on the orbit and the 
RD, as they may redistribute angular momentum among the different components. 



3.1. Torques acting on the NS 



The torques acting on the NS are the following. 



1. The spin-down torque due to the magneto-dipole losses (lLipunovlll992l ) 



N 



md 



2 Hi 



sin 2 x fei 



being k\ the direction of the NS spin, x is t ne angle between the magnetic and the 
spin axes, and /ii = B 1 R\/2 is the NS magnetic moment's modulus. In the following 
we assume k\ coincident with the unit vector k perpendicular to the orbital plane, and 
X = 7r/2, i.e. the NS's magnetic moment fi 1 lies in the orbital plane. 

We may assign a characteristic time scale to the torque (|BD 
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md | 



h 



10 45 gcm 2 y Vl0 33 Gcm 3 



f'i 



-2 



Pi 



10 ms 



(9) 



where P\ and I\ are NS spin period and moment of inertia. On account of the strong 
dependence t m d oc Pf , this torque is relevant in the early NS life, when Pi is small. It 
is fundamental in braking the initially fast NS rotation, but it soon becomes negligible 
and does not affect the final equilibrium rotation period. In the present analysis we 
assume a constant magnetic field, since a possible field decay does not affect our results, 
at least qualitatively. 

2. The torque due to the interaction between the magnetic moments of the NS and of the 
secondary star. The torque acting on the NS due to this interaction is 



N dd = tx x x B 2 (l) 



(10) 



Jackson! Il975l ). where -B 2 (l) is the RD magnetic field in correspondence of 



(see e.g. 

the NS: this can be written in terms of the dwarf's magnetic moment /x 2 as 

B /jn = 3 (/z 2 • n) n - /x 2 



(11) 



(|Jacksonlll975l ). where n is the unit vector pointing from the RD to the NS. We fix a 
non-rotating frame of reference centred on the NS, whose xy plane coincides with the 
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orbital plane. The torque ffTUl) also works on the orbital angular momentum to make 
n 2 an d ™ co-planar (IKing et al.lll990l ). We therefore assume that these vectors all 
lie in the orbital plane xy, so 



n 



cosi? 
sin •& 




COS01 

sin 0i 




COS 02 

sin 02 




(12) 



being d the orbital true anomaly, 0i and 02 the inclinations of the NS's and RD's 
magnetic moments with respect to the x-axis. After some algebra, the torque (1101) 
reads 



N 



dd 



a" 



2 cos(0 2 - 0) sin(0! - i?) + sin(0 2 - 0) cos(0! - 0) fc. (13) 



As a final remark concerning iVdd, we observe that it does not contribute to lock 
the period of the NS to the orbital value, but rather its phase angle 0i (see also the 
discussion in Appendix |AJ. Phase locking is a stronger condition than simple period 
locking. In the first case, the spin period and the orbital period coincide, but 0i may 
have any value. In the second case the orientation of 0i is fixed by the dynamics. 

Omitting the geometrical factor, we define the characteristic time scale of the dipole- 
dipole torque (fl3"j) 



tdd 



\N, 



7 x 10 fc 



<id 



10 45 gcm 2 / \10ms 



Pi 



10 67 cm f 



-i 



P. 



orb 



6.67 hr 



yr 
(14) 

This torque is weak when the NS rotates fast, but when its period is in the range of 
the hour, the time scale ( fl4|) drops to few years, and is therefore important for the 
final equilibrium. 

3. A dissipative torque acting to synchronise the NS spin with the orbital period. The 
nature and strength of this torque have been investigated by a hos t of autho r s, but 
there is still no general consensus about it. In their pioneering paper. Ijoss et al.l (119791 ) 
considered that the asynchronous orbital revolution of the secondary in the primary's 
magnetosphere would induce strong electrical currents in the secondary's atmosphere, 
which are dissipated by Ohmic decay, leading t he compact star t o rotate synchron ously 
with the orbit. In a series of papers Campbell (]Campbelllll983l . 11984 Il999l . 120051 ) con- 
sidered the role of turbulence in the secondary star, and concluded that the turbulent 
dissipation wo uld be much more effective than Ohmic decay in dissipating the cur- 
rents (but see iLambl Il985l ). Yet this torque becomes quite inefficient if ui\ > fi or b 
(fiorb = 2 7r/P or b is the orbital frequency), and therefore the compact star's initial spin 
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period must not be too different from the orbital period to reach synchronism. Other 
authors suggest different mechanisms for the dissipat ive torque such as unipolar in 



duction (jChanmugam and Dulkl Il983l : IWu et al.l 120021 ) or MHD torques ( iLamb et al. 
1983 ; Kaburaki 19861 ) . Despite their great differences, essentially all the models write 
the dissipative torque as 

AT diss = / ^, (15) 



with a typical time scale 



t. 



diss 



diss 



(16) 



necessary to synchronise the primary; here / ~ 1 is a model dependent dimensionless 
factor. The estimates of tdiss given by the diff erent models disagree by orders of mag 



nitude: typical values range from about 10 3 s (IKabu raki 



1986 



Chanmugam and Dulk 



19831 : IWu et al-lliooi ). to more than 10 6 yr JCampbelll ll983T) R In view of the large 
uncertainty about iVdiss, we have chosen to parametrise it. We would like to single out 
the dependence of this torque on the asynchronism between the NS's rotation and the 
orbital period, so we write 



. r T k>l — ^orb f 
iVdiss = ~ h K, 



(17) 



'syn 



where the synchronisation time scale r syn is a free parameter of our model. We will 
calculate the evolution of the NS's period for several values of r syn . An important 
caveat is in order here. The time tdiss strongly depends on the orbital separation a (see 
Equation (fl5l)). so the dissipation torque is inefficient in wide binaries. Essentially all 
the models for Polar CVs call for the condition Ra > a for synchronism. As the strong 
dependency on a is not explicit in the parameter r syn , we must check for consistency 
that in our synchronous solutions the condition Ra > a is fulfilled. 

If the dissipation torque is inefficient (i.e. r syn exceeds the age of RCW 103), the 
magnetic torques are unable to synchronise the NS, contrary to what happens in a 
Polar CV. In this case, IE would be more similar to a IP system, governed by material 
torques, rather than a Polar-like system ruled by purely magnetic interactions. 

4. The accretion and propeller torques. In its youth the NS spins too fast to a ccrete 
matter: this is the so-called "propeller" regime (see e.g. iRomanova et al.l 120031 ) . The 
mass coming from the companion star flies around the NS's magnet osphere, extracts 
angular momentum and spins the NS down. The exact form of the torque in the 



1 These times are meant as scaled from a white dwarf to a magnet ar. 
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prope ller stage is not well understood: we adopt the formula suggested by iLipunov 



(119921 . p. 159) 



N p = -kf§-k (18) 

-^co 

where kt is a constant to be determined later. 

This torque is very efficient, as shown by its short characteristic time scale 

t = ^~3xl0 2 AT 1 ( h \ ( r 2 (^-\ (^-) s 

p \N P \ 1 Vl0 45 gcm 2 ; U0 33 GcnrV \IAM Q ) VlOms; ' 

(19) 

The accretion torque depends on the mass accretion rate Mi on the primary, as well 
as on the arm's length of the material stress. In the case of a Polar, the flow from the 
secondary is directly accreted as soon as it crosses the inner lag rangian point Li : if the 



orbit is only moderately eccentric, we may adopt the formula (IPlavec and Kratochvil 



19641 ) 



b 1 = a (0.500 - 0.227 log 10 q) ~ 9.5 x 10 10 cm (20) 

to approximate the distance b\ between the inner Lagrangian point L\ and the NS 
centre. The torque's arm b\ yields 

N a = Mi (G M h) 1 / 2 k, (21) 

whose associated characteristic time scale is (up to factors of order unity) 

ta ~ L7 x 102 (lo^gcmO (lo^gs-i) {eWk) (e^Thr") yr ' (22) 

We describe both the accretion and the propeller regimes with a unique accretion/propeller 
torque 

AT ap = N p + N a . (23) 

We write Mi as a function of the mass transfer rate M? from the secondary. We 
postulate the relation 

M l = a\M 2 \- (24) 

where a ~ if N a <S N p and a ~ 1 if N a 3> N p . Since the accretion torque is dominant 
if b\ ^ R co and the opposite inequality is true if the propeller is dominant, we choose 

-wh' (25) 
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which may be rewritten 



a 



x Ma 



orb J 



where 



(26) 
(27) 

We determine the constant k t so that iV ap = when R co = b%, i.e. when the centrifugal 
barrier opens. After some algebra, we retrieve the formula 

1 A / U)\ 



\ = \(q) = (b 1 /a) 3 (1 + q) 
is a function of the mass ratio q only (cf. Equation (120jl ) 



N ap =\M 2 \ (GM60 1 / 2 



a (wi/n. 



orb J 



A 
2 



orb 



(28) 



In the accretion limit u;i <C A _1 / 2 fi or b the propeller torque is small, and this expression 
reduces to that of pure accretion (l2"Tj) . On the other hand, in the propeller regime 
ujx 3> A _1//2 f2 or b the first term in square brackets is negligible, and the torque reduces 
to pure propeller iV a p oc — uf, as in Equation ([TBI . 

Incidentally, we note that the NS spin equilibrium period P cq = X~^ 2 P rb defined by 
A ap = is the Keplerian period at L\ of a mass orbiting the NS, and c oincides with 
that predicted by the diamagnetic blobs accretion model put forward by King! (119931 ) 
and IWynn and King! (119951 ) for highly magnetised Cataclysmic Variables. 



The torques analysed in this Section are the building blocs for the overall torque acting 
on the NS: this will be made explicit in Section HI 



3.2. Torques Acting on the Secondary Star 



The reason why we consider also the torques on the secondary star is that they may 
couple to the orbit to alter P rb, thus affecting synchronism (see also lKing et al.lll990l for a 
fuller discussion). In general, the torques on the secondary star are: 



1. The magnetic braking torque. If the RD is magnetic and blows a wind, this is coupled 
to the magnetic field and may extract angular momentum from the secondary. This 
torque only depends on the RD period and not on the orbital one, so magnetic braking 
does not directly affect the orbital period. However, it may couple with the tidal torque 
(discussed immediately below) to transfer angular momentum from the orbit to the 
RD. 
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2. The tidal torque reduces the relative rotation of the RD and the orbit, and of course 
vanishes if the secondary rotates synchronously. This torque may also redistribute 
angular momentum between the orbit and the RD. Stars with convective envelopes 
(which is most likely to be the case of the our companion star) in binary systems with 
orbital periods of few hours have tidal synchronisation times ~ 10 2 yr (Equation [T]), 
so we shall assume that the RD is always locked to the orbital period. The magnetic 
braking extracts angular momentum from the secondary, which is then taken away 
from synchronism. The tidal torque, on the other hand, is very effective in locking 
the secondary back to orbital synchronism. As a consequence, the synergic action of 
these torques pumps a ngular momentum from the orbit to the secondary, keeping it in 



synchronous rotation. (jVerbunt and Zwaanlll98ll ) 



If the secondary is not synchronous, there are two effects. First, the dipole-dipole 
torque (fT3l keeps oscillating from positive to negative, so it cannot lock the NS's phase. 
Second, the secondary can exchange some if its intrinsic angular momentum with the 
orbit. In view of our estimates (Section 13.31 below), we do expect that the amount 
of transferred angular momentum is too small to affect the orbit. Summarising, the 
secondary's possible lack of synchronism does not directly alter the primary's period 
locking, but may invalidate its phase locking. 

3. The dipole-dipole torque. This is exactly the same kind of torque ( |T3|) acting on the 
primary: its expression is readily derived by swapping 0! and 02 in Equation (1131) : 



iV dd — 3 



2 cos(0i - ■&) sin(0 2 -■&) + sin(0! - 0) cos(0 2 - ■&) k (29) 



This torque depends on the true anomaly as well as on the stars' magnetic moments, 
and is able to exchange angular momentum between the secondary and the orbit. The 
synchronisation time scale for this torque is 

a = ^A^ 3xl0 9 ( h \ ( p * V 1 ( ^ r 1 ( Porb V vr 

l-W'ddl V10 54 gcmV V6-67hr ^ Vl0 67 G 2 cm 6 / \Q.Q7Yir) y ' 

(30) 

where P 2 is the period of secondary. Since t' dd is much longer than the age of RCW 103, 
the RD phase is not likely to be locked, although its period is most likely synchronous 
with the orbit. 

In summary, we assume that the secondary's rotation is synchronous with the orbit. 
The instantaneous phase 02 (i-e. the position of the RD magnetic moment /^ 2 ) is 

02 = ^orb t + 02, (31) 

where the phase 2 is arbitrary due to the lack of phase locking. 
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3.3. Orbital Torques 

In this Section we derive and discuss the torques acting on the orbit of the binary system. 
Since some torques and the mass transfer alter the orbital period, they must be taken into 
account to study the synchronisation. As usual, we suppose that the orbit is circular, so the 
orbital angular momentum reads 



Go 



1/2 



Jorb = I ~w) Ml M2 k (32) 

The rate of change of J orrj is due to the action of all the stellar torques containing the true 
anomaly $ (or the orbital frequency fi or b), plus the tidal/magnetic braking torque N mag due 
to the angular momentum exchange with the secondary discussed in Section 13.21 : 



orb 



-N dd - N' dd - N ap - N diss - iV mag . (33) 



Gravitational radiation has not been included: it is generally important for systems with 
short period (P Q rb ^ 2 hr), and it is immaterial here. 



The torque iV ma g is (jVerbunt and Zwaanlll98ll ) 



iV mag ~ 0.5 x 10" 28 I 2 R\ Q 3 orh (34) 
(in cgs units). This torque drags orbital angular momentum on the time scale 

t mag = J orb /iV mag * 10 8 (-^-] (-^-) 10/3 (-rJ^-X 1 W, 



0.3 M J V6-67hry V 1054 g cm2 / V3-3xl0 10 cm 

(35) 

where m is the reduced mass of the binary system. Since this time scale is very long, the 
orbital angular momentum yield to the secondary may be neglected. 

Ignoring the exchanges of orbital angular momentum with the RD, we focus on the 
exchanges with the rapidly spinning NS. Adopting the parametrisation ( 1241) for the mass 
accretion rate, the logarithmic derivative of the modulus of Equation (1321) gives 

Jorb = i a _ jMal 2(1- qq 2 ) + q(l - a) 

J mh ~~ 2 a M 2 2(l + q) 1 j 

Plugging the expressions of the relevant torques into ( 1331) . and using the Kepler relation 

-2/; 

orb 



2/3 

a oc fi orb we retrieve our equation for the evolution of the orbital frequency 



^orb 3 x% 

^2 ma 2 O, 



orb 

\M 2 \ (GMbj) 1 / 2 
h 



sin(0 1 + 2 -2^) + ^ 

'syn 

l h /,, / n ^2 

l+A KM, rb )2 

\M 2 \ 2(l-qg 2 
2 M 2 1+q 



— |A (<^i/0 or b) 



+ 



3 \M2\ 2(l-gg 2 )+g(l-g) ( o 7 \ 
2 Mo 1+a \°'J 
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We briefly discuss the magnitude of the terms on the right-hand side of Equation (j3*7|) . 
The terms in square brackets are the same dipole-dipole, dissipative and accretion/propeller 
torques (with allowance of signs and geometrical factors) acting on the NS; all these terms 
are multiplied by the small factor ( = 3 Ii/ma 2 Q or \, ~ 10 -6 , showing that the time scales 
for the orbital evolution are longer than the time scales for the NS' spin evolution by a factor 

C _1 ~ io 6 . 

Mass transfer unbalances the relative masses of the two components, and only becomes 
significant over the long time scale 

h * M 2 /\M 2 \ = 2.5 x 10" (J^j (jol^L) (38) 

In summary, the orbital torques are characterised by time scales much longer than the 
current age of RCW 103; the orbit may be then be regarded as fixed, and in the following we 
only consider the torques on the NS. The results presented in the next Section are based on 
this hypothesis; however we had an additional check of its validity by running a numerical 
calculation where the orbital parameters were allowed to evolve according to Equation (|37|) . 



4. Evolution of the NS Rotation 



We adopt (the phase of the NS magnetic moment) as the NS Euler free rotation angle 



(e.g lLandau and Lifshitzlll969l ). Putting all the torques (IE]), ( TTBl . ( fTTl) and ( |28l) together, the 
equation of motion for the NS's phase angle 0i reads 



n3 



2 cos(02 — $) sin(0i — $) + sin(0 2 — ■#) cos(0i — $) 



2 l4 3 



T syn 



l+A (a;i/r2 orb )2 



| (^l/Oorb)' 



(39) 



where the true anomaly = fi or b t is assumed to be zero at t = 0. As noticed above, we keep 
the orbital parameters constant, as their variation due to the angular momentum exchanged 
with the stars is negligible. 

As shown in Appendix [A], Equation ( |39i) admits an asymptotic constant solution, i.e. 
a stable NS equilibrium spin period. In the absence of material torques (M2 = 0), this 
period is attained thanks to the second and the third term in Equation (|3"9"]1 . leading to full 
synchronism with the orbit. The magneto-dipole torque brakes the initial fast NS rotation, 
and later on the dissipative torque brings the NS spin to orbital synchronism. In presence 
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of material torques, it is the interplay between propeller and accretion torques that drives 
the young magnetar into equilibrium (as implied by Equation ([6])). In this case, the orbital 
period does not coincide with the NS's spin equilibrium P eq , but it is quite close to it 
(P eq ~ 0.5 — 0.7 -Port.)- Even in this case the NS magnetosphere keeps filling the binary 
Roche lobe, warranting the consistency of our results. 

We integrated Equations (l39l) with a step-adaptive fourth order Rosenbrock algorithm 



( iPress et al.ll2002l ). This choice is motivated by the possible stiffness of the problem due 
the diverse time scales of the torques. The NS and RD magnetic moments are the free 
parameters of the model, together with \M 2 \ and r syn . The magnetic moments of the NS 
and the RD are \i\ = 10 33 G cm 3 and fi 2 = 10 34 G cm 3 respectively. The NS's moment of 
inertia has the standard value of I± = 10 45 g cm 2 . We let the mass transfer rate and the 
synchronisation time to vary in each run, in order to explore different evolution paths. The 
initial phases <fi\(t = 0) and 02 (t = 0) have been randomly chosen, with an initial NS rotation 
period of 3 ms. The RD has been always considered synchronous. In all our runs the NS 
final equilibrium spin period has been set equal to 6.67 hr (see Appendix [X]), and we have 
derived the orbital period in a self-consistent manner. 

For consistency, we have also checked that the condition that Ra > a, is always met in 
our runs. The NS Roche lobe is dominated by the star's magnetic field, and the mass from 
the secondary is effectively channelled from L\ to the star's polar caps without forming an 
accretion disc, consistent with the torque (l28i) used in our calculations. 

It is convenient to present first our result in the case where the mass transfer rate 
M 2 vanishes. The material (i.e. accretion and propeller) torques are always zero, and the 
evolution is controlled uniquely by the magnetic torques. Figured] plots the evolution of the 
NS' period for different values of the synchronisation time: r syn = 10 yr, r syn = 100 yr and 
r syn = 1000 yr, respectively. 

The dissipative magnetic torque plays, in this key role in driving the NS spin 

into orbital synchronisation, attained within ~ 2000 yr if r syn is not longer than few hundred 
years. A magnetar-like field and a close orbit are necessary and this i s cons i stent with 



estimates given in the literature, possibly with the exception of ICampbelll (119831 . 120051 ) (see 
also the discussion in Section 13.11) . 

The situation changes if the mass transfer rate is non-negligible. Figure [2] plots the 
evolution of Pi and of the synodic normalised frequency 

£ = (ux - Q OTh )/Q orh , (40) 

for |M 2 | = 10 13 g s -1 . Here the material torques are relevant, and lead to the equilibrium spin 
period within few hundred years, fairly independent on the exact value of r syn . The orbital 
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period is longer than the NS's equilibrium spin period; we find P or b = 6.9 hr if r syn = 10 yr, 
-Porb = 8.7 hr if r syn = 100 yr and P or b = 11.5 hr for r syn = 1000 yr. An additional run with 
r syn = 10 7 yr shows that the magnetar attains its spin equilibrium within ~ 10 3 yr, with an 
orbital period of ~ 12 hr, so that in all cases we find P sp m/ -Porb ^ 0.5. In magnet ic CVs period 



ratio s larger than 0.5 correspond to stream-fed accretion onto the white dwarf ([Norton et al. 



20041 ). Consistently, we find also here that the condition Ra > b\ is always met: the NS's 



Roche lobe is filled by the magnetosphere, and no accretion disc can form. In addition, 
for T syn < 1000 yr, the stronger condition Ra > a is also met: this warrants that this 
(parametric) synchronisation time is consistent with the requirement of a strong magnetic 
interaction between the stars (Section l3.ip . 

In the same Figure [2] we show the evolution of Pi and £ for the much larger mass 
transfer rate IM2I = 10 15 g s _1 , corresponding to the peak X-ray luminosity observed in IE. 
This value should provide an upper limit to the influence of the mass accretion rate on 
the period evolution. For this value of \M 2 \, the condition Ra > b\ is violated, but since 
Ra is still larger than the circularisation radius P C irc of the accreting fl ow, the transferred 



mass is unable to circularise and form a disc (see e.g. iFrank et al.ll2002l ). For such a strong 
material torque, the equilibrium period is attained in less than ~ 10 yr during the early 
"propeller" phase which extracts a great amount of angular momentum from the spinning 
young magnetar. 



5. Discussion 

In the previous Section we have shown that if IE hosts a highly magnetic NS orbiting 
a RD, the NS is spun down rapidly to an asymptotic equilibrium period, either by magnetic 
and/or material torques. The magnetar-like NS loses memory of the spin acquired at birth, 
and approaches a spin period close to the binary orbital period. 

If the material/propeller torques have been small for an appreciable fraction of the 
system's age, then magnetic torques have been dominating since birth, and orbital synchro- 
nism is attained within 2000 yr if r syn < 100 yr. The X-ray emission in this case should be 
dominated by mechanisms similar to those seen in Magnetars (Section [2} . 

On the other hand, if the material torques have been important, IE is more similar 
to a young accreting low-mass binary. The close resemblance of the light curve of IE with 
Polars and Intermediate Polar Cataclysmic Variables suggests that mass transfer currently 
takes place in IE via filling of the RD's Roche lobe. (Note that this does not exclude the 
action of magnetic torques in the earlier life of the source). Although the magnetic torques 
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are unable to synchronise the system, the hypothesis of a magnetar-like field is still essential 
here, since the NS magnetosphere must dominate the star's Roche lobe. If B\ ~ 10 15 G and 
the orbital period is of 6 — 12 hr, then Ra > b\ and the flow from the RD provides a material 
torque with arm of length ~ b±, which brings the NS to the observed equilibrium spin period 
P eq = 6.67 hr within ~ 100 yr, independent of the magnetic torques. The orbital period 
in this case still depends on r syn , and is longer than P e a- m the limit of a very inefficient 
dissipative torque we find P or b ~ 2P eq (i.e. P or b — 12 hr) q 

A remark about our use of the circular approximation. As can be seen from Equa- 
tion which rules the secular evolution of the NS's spin, the circular approximation 
enters: i) the distance between the components in the dipole-dipole interaction ; ii) the 
orbital frequency f2 or b in the dissipation term and the accretion/propeller torque ; iii) the 
arm b± (the distance between the neutron star and the inner Lagrangian point), which is 
calculated in the Roche formalism. Note that (whatever the eccentricity) point iii) is imma- 
terial if the RD fits within its Roche lobe. In this case, there is no mass transfer, and the 
accretion/propeller torque vanishes; if the companion is large enough to allow mass accre- 
tion, the arm h\ only enters as a weak power in the torque, and a possible deviation from the 
circular geometry is unlikely to affect this torque to a significant amount. As far as points i) 
and ii) are concerned, the introduction of a non-zero eccentricity would induce an orbital 
modulation in the dipole-dipole and in the dissipation torque. This modulation is not able 
to alter the qualitative results presented in this paper. 

In this model, IE would be the first young low mass binary observed in X-rays in the 
first ~ 2000 yr of its evolution, and there might be a connection between the magnetar 
nature of the NS and this premature onset of an X-ray phase if it comes from accretion. 
Two questions arise: 



2 It is worth spending here few words on the upper limits of the orbital period. If the companion star 
tightly follows the mass-radius relation valid for low main sequence stars, then the orbital period is unlikely 
to exceed by much P or b ~ 12 hr, since this binary separation brings the NS to its spin equilibrium period 
due to the sole action of the material torques. These imply a Roche-lobe filling donor star, which in our 
calculations we assumed to be of 0.4 M©, close to the upper limit set by the optical/IR observations. Were 
the system wider than this, the star should fall within its Roche lobe, making mass transfer impossible, as 
well as the action of the material torques. On the other hand, mass transfer is possible if the companion 
star is larger than a main sequence star of the same size, e.g. because it absorbed energy from the early 
relativistic wind emitted by the NS. In this case the limits on the orbit's size are more uncertain, but by 
analogy with the period distribution of bona fide Polars and IPs, we may estimate that the orbital period 
is unlikely to exceed by much P or b ~ 12 hr. The current observations, on the other hand, are much less 
constraining. Our X-ray data do not yield any evidence for a periodical modulation larger than 6.67 hr. In 
any case, it would be quite difficult to detect a modulation at a few percent level and a long periodicity 
(10 — 15 hr) superimposed on the very strong pulsation at 6.67 hr. 
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1. Does the pulsar wind of the young magnetar induce a mass loss from the RD? Can this 
wind vaporise the RD? How does this effect constrains the initial secondary's mass, or 
the NS initial period? 

2. Looking ahead, what is the subsequent evolution of IE? 

As far as the first question is concerned, the strong wind of photons and relativistic 
particles of the magnetar may perturb the RD thermal equilibrium, and may in principle 
cause its partial evaporation. The problem we address here differs from the evaporation of a 



main s e quence star by a millisecond pulsar in a narrow binary system (e.g. iRuderman et al. 



1989aB l Tavanilll992l ). In these systems, the NS magnetic field is low (10 8 — 10 9 G) and its 
weak irradiation lasts for more than 10 8 yr. In a millisecond pulsar environment, the main 
sequence star has enough time to adjust to thermal equilibrium during the evaporation. 
In the present case this assumption is most likely invalid, since the irradiation from the 
magnetar is much more intense and short-lived (as shown below). 

If the magnetar in IE has a birth period P = 2tt/ujq of few milliseconds, it holds 
a rotational energy E TOt ~ 2 x 10 51 (Po/3 ms)~ 2 erg, which is dissipated mainly through 
magnetic dipole losses in the early years, when the propeller torques are still ineffective. The 
RD has a binding energy E RD ~ GM$/R 2 ~ 10 48 (M 2 /0.4 M ) 2 erg, so the RD needs to 
absorb only ~ 0.1%E rot to be evaporated. 

We can estimate the mass loss rate from the RD by comparing the NS energy loss (and 
impinging on the RD) with the kinetic energy necessary to sustain a wind, regardless of the 
physical interaction process [E In more detail, the spin-down luminosity of the magnatar 
decays with time as 

w*) = 5^ fi + — («) 

2 V Tmd/ 



3 The wind of particles fro m an ordinary milliseco nd pulsar is composed by energetic photons (X-rays 
or 7-rays) and TeV e^ 1 pairs (jRuderman et al.lll989af ). In the case of a magnetar, the occurrence of high 
energ y pairs is debatable, since the strong magnetic field may prevent their formation due to photon split- 
ting ( Duncan and Thompson 1992h . With this caveat in mind, we estimate how deep the impinging wind 
penetrates the companion star. Photons are stopped in the stellar atmosphere, having a mean free path of 
few cm in a gas with density ~ 1 g cm -3 . The electrons impinging on the companion star yield their energy 
chiefly for Inverse Compton losses on the radiation field, and for relativistic Bremsstrahlung. The Inverse 
Compton and Bremsstrahlung stopping lengths are lie — 10 6 cm and Zbs — 10 5 cm, respectively (|LansJll98Q 
p. 467): since Zbs and he are both small in comparison to the stellar radius, essentially all the energy carried 
by the relativistic wind is deposited in the star's envelope. 
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where r md is the magneto-dipole decay time 

= - 4 -6 x 10 3 s (jj^) " (^-) 2 • (42) 

4 h(uJq V lO 1 " G cur / \ 3 ms / 

Equations (|4*T|) and (14*2"]) show that (i) the magnetar is braked to a period of 1 sec after ~ 16 yr; 
(ii) most of the rotational energy is dissipated in the first year; (iii) the spin-down luminosity 
decays by a factor ~ 10 10 over this very short time. At P ~ 10 sec the magnetar in IE brakes 
under the action of the magnetic and material torques, but at this period the bulk of its 
rotational energy has already been dissipated, and is but a tiny fraction of the companion's 
binding energy. For this reason, any later mechanism that brakes the NS and which may 
- in principle - transfer energy from the NS to the companion, is most likely to have little 
effect on the latter. 

We suppose that this relativistic wind is emitted isotropically, so a fraction (R 2 /2a) 2 
hits the companion star; a fraction r\ of this impinging power ext racts an evaporation wind 



from the RD. Following Ivan den Heuvel and van Paradijsl (119881 ). we equate the effective 



power of the relativistic wind with the evaporation wind to estimated M 2 . 

, , x9 G M 2 dM 2 , , 

77 (R 2 /2a) 2 L SD = — ^ (43) 

where we have assumed that the evaporation wind speed is close to the escape velocity from 
the star's atmosphere. We obtain the differential equation 

GM 2 dM 2 IfRA 2 hj4 ( 1 + J-Y 2 (44 ) 



R 2 dt 2 \ 2a / T md V T md 

for the evolution of the star's mass M 2 (t), which can be solved once we have specified the 
relation between M 2 and the star's radius R 2 . This relation is uncertain, so we consider two 
limiting cases aiming to bracket the actual behaviour. 

• On the one hand, we assume that R 2 is given, instant by instant, by the relation (|7j) 
valid for lower main sequence stars. This is the most optimistic limit, since it assumes 
that the companion star has enough time to reduce its size as a reaction to the evapo- 
rating flux. This may be not very realistic, however. The star reacts to the relativistic 
wind on the Kelvin-Helmholtz time for the envelope involved in the energy deposition 

tkh^ttt^I^IOV, (45) 
it 2 1j 2 rC 2 



which is much longer than the magneto-dipole time scale (142]) ; L 2 is the RD's luminos- 
ity, and Ibs is the penetration depth of the relativistic wind, where thermal equilibrium 
is perturbed (see note). 
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• On the other hand, we may suppose that the star's envelope swells up, filling its own 
Roche lobe. In this case, we assume for R2 the size of the Roche lobe at the evaporation 
onset. This is the most pessimistic limit, since in this case the section presented by 
the star to the impinging flux is the largest, resulting in a fast evaporation rate. 

Figure [H] shows the final mass of the companion star versus its initial mass in these two 
limiting cases for several values of the parameter 77. The fixed parameters are the NS's initial 
period (3 ms) and the orbital separation (a = 1.5 x 10 11 cm), although an increase of a is 
expected on account of the mass lost by the system. 

In both cases, the secondary's final mass is M 2 < 0.4 M provided that 77 < KT 1 , for 
an initial mass M2 < 0.6 M . For efficiencies 77 > 10 _1 and initial masses M2 < 0.6 M Q the 
secondary is completely evaporated in the Roche lobe filling case. If the M2 — R2 relation (J7J) 
is followed, M 2 drops below 0.1 M & if 77 > 0.1. We may conclude that (in either scenario) 
the secondary is left with a mass M 2 < 0.4 M & if its initial mass is M 2 — 0.6 M Q and the 
efficiency is rj ~ 1 0~ 2 ; this valu e is consistent with the estimates rj ~ 10~ 2 — 10 _1 cited in 
the literature (e.g. 
by about ~ 10%. ( 



e.g. 



Tavanilll992l). The mass lost dur ing the evaporation widened the orbit 



Postnov and Yungelsonl 120061 1 



The evaporation scenario we have explored is perhaps the most unfavourable for the 
survival of the companion, as some assumptions may be a little bit unnecessarily pes- 
simistic. First of all, we have assumed that the relativistic wind is isotropic. This may 
not be the case: in some models for the birth of magnetars (in connection with the problem 
of lo ng Gamma-ray Bursts) most of the rotation energy is removed by highly collimated 
jets (jBucciantini et al.ll2007l ). If the collimated beam keeps away from the RD, the problem 
of evaporation disappears. Second, we assumed a fast initial NS s pin period {Pn = 3 ms), 
as required to build up the magne tar's field by a a - dynamo (jDuncan and Thompson 



1992 



Thompson and Duncan! Il993l ). There are, however, alternative models which explain 



the formation of strong magnetic fields as du e to the amplification of high fossil fields due 
to flux conservation during the core collapse (jFerrario and Wickramasinghd 120051 . 120061 ); in 
this case there is no need of a very NS fast initial rotation period. On account of the strong 
dependency of Lsd oc Pq 2 (Equation ( 14TI) ). also the wind hitting the RD is less energetic. To 
summarise the discussion relative to point (1), there are several possible scenarios in which 
the evaporation of the companion is not an issue. Even in the worst case we have analysed in 
some detail, for reasonable combinations of the parameters the RD may survive evaporation. 

Although the early relativistic wind is unable to evaporate the whole star, it is not with- 
out consequences. Indeed, the secondary star quickly recovers its hydrostatic equilibrium, 
but takes a longer time (see Equation fT45]) ) to restore its thermal equilibrium. It is quite 
difficult to address the structure of the companion after it has been shaken by the magnetar 
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wind. The companion is likely to be larger than a ZAMS star with the same mass, so it may 
fill its Roche lobe and transfer mass to the NS, powering (part of) its X-ray luminosity. 



As far as query (2) is concerned, we expect that when the RD returns to thermal equilib- 
rium, it shrinks, fitting inside its own Roche lobe. Mass transfer ceases and only the residual 
magnetar like emission remains. This likely occurs on the Kelvin-Helmholtz time (I45p . Later 
on, we expect that also the magnetar field will decay. Current theoretical models predict that 
the present field of 10 15 G will decay to IP 12 G (typical of the magnetic field of NSs in ordinary 



LMXBs) in < 10 8 yr flHevl and Kulkarnilll998l : ICofpi et al.l 120001 : iGeppert and Rheinhardt 



20021 ) on account of non-linear field effects like ambipolar diffusion and Hall drift. Later 
on, accretion f rom the companion sets in that can induce further field decay, driy ing it be- 
low ~ 10 10 G ( IBhattacharya and van den Heuvell Il99ll ; IGeppert and Urpinl Il994l ). As the 
companion star will exhaust its nuclear fuel in few 10 9 yr, it will swell up as a giant, filling 
its Roche lobe and transferring mass to the NS. Since the NS's field is weak by now, the 
accretion flow may settle in a disc. The ensuing spin up torque may accelerate the NS up to 
few milliseconds, and at the end of this stage the system would look like a bona fide LMXB. 

A speculative deduction of our model is that at least some LMXBs may have systems 
like IE as progenitors. The magnetic field of a NS born as a magnetar in a low mass binary 
might keep stronger than the field of a NS in a LMXB or of a binary millisecond pulsar. This 
kind of higher field is observed in the seven Soft X-Ray Transients LMXBs where an X-ray 
pulsation is detected during type I burst. This will be explored in a forthcoming paper. 



6. Summary 

In this paper we have put forward a model for the unique X-ray source 1E161348-5055 
(IE) an the centre of the young supernova remnant RCW 103. In our model IE is a binary 
system made by a NS and a lower main sequence star. We interpret the signature of 6.67 hr 
as the NS's spin period. This requires that the NS is endowed with a strong magnetic field 
of ~ 10 15 G. 

The spin evolution of the NS depends on the relative importance of the material and 
magnetic torques. If the magnetic torque is dominant (which is the case if most of the X-ray 
luminosity is powered by magnetar emission), the NS's rotation is almost synchronous with 
the orbital period, similarly to what happens in Polar Cataclysmic Variables. On the other 
hand, if the material torques are more relevant, then the spin period is about 50 — 70% of 
the orbital period. In this case, the system would be similar to an Intermediate Polar CV, 
with a NS instead of a white dwarf. The magnet ar-powered emission and the full orbital 
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synchronisation scenario is most likely to o ccur if the secondary s tar is small (i.e. significantly 



below the current upper limits derived by iDe Luca et al.l (120061 )). 



Alternatively, an early relativistic wind from the NS may have stirred an evaporation 
wind from the companion. On account of the large amount of energy injected into the RD's 
envelope, the RD may be out of thermal equilibrium. It may be larger than a main sequence 
star with the same mass, possibly powering by accretion the NS' X-ray luminosity. This 
phase is expected to last 10 4 yr, the Kelvin-Helmholtz time scale of the RD layers involved 
in the energy deposition from the early relativistic wind. As soon as the star recovers its 
thermal equilibrium, the star fits again within its own Roche lobe, and mass transfer will 
end. 

If the model we have suggested is correct, then IE is the first case (to our knowledge) of a 
magnetar hosted in a low mass binary system. This may spur some interesting developments 
in the studies of the evolution of narrow binary systems harbouring a NS. 
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A. Equilibrium Solutions 

In this Appendix we study in more detail the properties of the equilibrium solutions 
found numerically in Section HI To this end it is convenient to cast Equation ([3U]) in a 
dimensionless form. We measure the time in units of f2~ b , i.e. we set t = t £l or b', we then 
introduce the synodic phases 

0'l = 01 " t, 2 = <P2- I (Al) 

If the secondary is synchronous, then 2 is independent of time. After some algebra, Equa- 
tion ( 1391) reads 

§ = e W(£) (A2) 
at 

where 

e = A (A3) 

£ = | £ 2 + n dd (sin 2 sin^ — 2 cos0 2 cosc/4) (A4) 

W(0 = -n md (1 + 3 - n d]ss £ + n ap [ 1+A( \ +g)2 - § (1 + £) 2 ] , (A5) 

and 

n md = f 4^ (A6) 

n d iss = 1/ (O orb r syn ) (A7) 

n ap = ^fg 61)1/2 (A8) 

^ orb 

n dd = (A9) 

are the dimensionless torques on the NS. Equation (1A21) has the form of an energy equation 
in a non-conservative system, on account of the forcing term £ W(£) on its right-hand side. 
The dipole-dipole torque does not appear in the braking function W, showing that this 
torque does not help to attain synchronism, since it does not alter the "energy" S. 

If there is no mass transfer, then n ap = 0, and since n md <C ridiss Equation (1A2jl reads 
approximately 

^ ~ -n diss £ 2 (A10) 
at 

Starting from a non synchronous configuration (£ ^ 0), the right-hand side of this equation 
is negative, so dS/dr < 0. The system then evolves towards the minimum of £ , given by 
£ = <p[ = and 

,/ 2 cos02 . ,/ sin 2 
cos 0i = — sin 0i = All 



v/1 + 3 cos 2 02 ^/l + 3 



cos 
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This is the synchronous solution, where the synodic phase <\>' x is locked to the angle defined 
by the last couple of equations. It is also interesting to note that the "potential" 

V(<f>i) = ridd (sin 02 sin^ — 2 cos0' 2 cos 01 ) (A12) 

in Equation ( 1A4I) is proportional to the the magnetic interaction potential energy between 
the stellar magnetic dipoles U = —fix ■ -62(1). The equilibrium synodic position angle of the 
synchronous solution therefore minimises the interaction magnetic potential energy between 
the two stars. 

The qualitative pattern of the evolution towards synchronism does not change if W 
remains negative. The only positive term in W which may change the sign of W and alter 
this conclusion is the first term in square brackets proportional to n ap . This is the accretion 
term, which tends to spin up the NS. We now show that for a large accretion rate, the 
NS tends to a stable asynchronous uniform rotation. If n ap exceeds a threshold value, the 
equation W(£) = has a root £o- Since W'(£) < for any £ > 0, this root is also unique. In 
addition, W(£) > if £ < £ and W(£) < if £ > £ - These results show that d8/dt < for 
£ > £0 an d dS/dt > for £ < £ : the system evolves towards the value of £ characterised by 
£ = £0 an d X = Co t + 0o- Since £ 7^ 0, this is an asynchronous solution: the system cannot 
attain synchronism when the accretion spin up torque exceeds a criti cal value, in agr eement 



with what is well known in the literature for Polar systems (see e.g. ICampbelllll986l ) 



The result is shown in Figure HI which plots the equilibrium synodic frequency £ against 
the mass transfer rate \M 2 \, for different values of the synchronisation time r syn . As expected, 
£ grows with \M 2 \- Below the critical threshold \M 2 \ ~ 10 5 gs _1 the equilibrium solution is 
always synchronous. As the dissipative torque decreases (i.e. for long r syn ), synchronism is 
more difficult to achieve. For \M 2 \ — 10 13 gs _1 (the mass accretion rate inferred from the 
source's X-ray luminosity) and the most favourable case r syn = 10 yr we find that the NS 
has a degree of asynchronism of a few percent. 
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Fig. 1. — Evolution of the NS's rotation period without mass transfer (M 2 = 0). The dotted 
line refers to the synchronisation time r syn = 10 yr, the dashed line to r syn = 100 yr and the 
dot-dashed line to r syn = 1000 yr. In all cases the orbital period is P or b = 6.67 hr. 
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Fig. 2. — Evolution of the NS's rotation period (upper panel) and of the normalised synodic 
frequency £ = {uj\ — VL or ^)/VL orh (lower panel). The two groups of lines refer to the mass 
transfer rates | Af 2 | = 10 13 g s _1 and \M 2 \ = 10 15 gs _1 . The lines' coding is the same as 
in Figure [TJ For |M 2 | = 10 13 gs~ 1 the orbital periods are P or b = 6.9 hr (for r syn = 10 yr), 
P orb = 8.7 hr (for r syn = 100 yr) and P orb = 11.5 hr (for r syn = 1000 yr). For |M 2 | = 10 15 g s" 1 
the orbital periods are P or b = H-4 hr (r syn = 10 yr), and P or b = 11.9 hr (both for r syn = 100 yr 
and r syn = 1000 yr). Note that for the mass transfer rate |M 2 | = 10 15 g s -1 the condition 
Ra > b\ is violated (see text for more details). 
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Fig. 3. — The final mass of the RD after the evaporation process from the mageto-dipole 
radiation from the NS versus the initial mass for different values of the efficiency rj. Upper 
panel: the radius of the secondary is always taken as the ZAMS radius for a lower main 
sequence star. Lower panel: the radius of the secondary is taken as the Roche radius at the 
evaporation onset. Dotted line: rj = 1CT 3 ; dashed line: rj = 1CT 2 , dot-dashed line: rj = 1CT 1 . 
In both cases, the orbital separation is fixed to 1.5 x 10 11 cm, correspondent to an orbital 
period of 6.67 hr for a binary system with M x = 1.4 M & and M 2 = 0.4 M Q . 
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Fig. 4. — The normalised synodic equilibrium frequency £ = (u>i — f2 or b)/fi or b plotted against 
the mass transfer rate |M 2 | from the secondary star, for different values of the synchronisation 
time r syn . Dotted line, r syn = 10 yr; dashed line r syn = 100 yr and dot-dashed line r syn = 
1000 yr. As |M 2 | increases, the value of £ tends asymptotically to the value £ ~ 0.7. 



